REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OM?No.  0704-0188 


PanM-wnrk  RcdiicdoTi  Proicct  ((T704-0188).  Washin^cc,  DC  20^3.  _ _ _  i\^  . 


Monageroent  Badgit.  Paperwork  Reduction  Project  (0704-0188),  Waahipgtcc,  E)C  20S0: 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

June  7,  1995 


3.  REPORT  TYPE  AND  DATES  COVERED  r::  _ 

August  1,  1992  to  Jandary  1! 


4.  TITLE  AND  SUBTITLE 

Direct  Numerical  Simulation  of  Acoustic-Flow 
Interactions  in  Solid  Rocket  Motors 


-  5.  FUNDING  NUMBERS 

F49620-92-J-0451 


6.  AUTHOR(S) 


S.  Mu,  S.  Revert,  and  S.  Mahalingam 


The  Regents  of  the  University  of  Colorado 
Campus  Box  19 
Boulder,  CO  80309-0019 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
Bolling  Air  Force  Base,  DC,  20332-6448 


AFOSR-TR-95 


05 


153-7118 


10.  SPONSORING  /  MONITORING  AGENCY 


1 2a.  DISTRIBUTION  /  AVAILABILI I  Y  b  I  ^ 


Approved  for  public  release, 
distribution  unlimited 


kW  d\d\ 


2b.  DISTRIBUTION  COD 


13.  ABSTRACT  (Maximum  200  words) 


A  summary  of  research  on  our  studies  of  acoustic-mean  flow  interactions  in  solid  rocket 
motors  is  presented  in  this  Final  Technical  Report.  Two-dimensional,  time-dependent 
solutions  of  the  compressible  Navier-Stokes  equations  have  been  obtained  in  a 
rectangular  duct  under  conditions  of  monochromatic  acoustic  forcing.  The  existence  of  a 
cut-off  frequency,  the  appearance  of  oblique  waves,  and  the  behavior  at  resonance 
compare  well  with  previously  published  analytical  predictions.  The  thickness  of  the 
acoustic  boundary  layer  and  its  response  to  imposed  disturbances  are  also  in  good 
agreement  with  theory.  The  code  is  being  modified  to  include  nonpremixed  combustion, 
which  will  allow  studying  acoustic  interactions  with  a  chemically  reactive  mean  flow. 


17.  SECURITY  CLASSinCATlON  OF  18.  SECURITY  C^SSIHCATION  OF  TfflS  19. 

report  page  I  r  ABSTRACT 

uhtLfutrer^  IJ _ 


15.  NUMBER  OF  PAGES 

41 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 

Standard  Form  298  (Rev.  2-89) 


19950821  056 


DTIC  QUALITY  IK'SPECTBD  B 


Direct  Numerical  Simulation  of  Acoustic-Flow  Interactions 


in  Solid  Rocket  Motors 


Siming  Mu,  Stephen  Hevert  ajid  Shankar  Maiialingam 


Final  Technical  Report 
AFOSR  Grant  No.  F49620-92-J-0451 
Contract  Monitor:  Dr.  Mitat  Birkan 

June  7,  1995 

CCR  Report  No.  95-02 


Center  for  Combustion  Research 
Department  of  Mechanical  Engineering 
University  of  Colorado  at  Boulder 
Boulder,  CO  80309-0427 


ABSTRACT 


A  summary  of  research  on  our  studies  of  acoustic-mean  flow  interactions  in  solid  rocket 
motors  is  presented  in  this  Final  Technical  Report.  The  objective  of  this  research  was  to  study 
the  effects  of  acoustic  refraction  on  the  mean  flow  with  imposed  acoustic  disturbances  of  vari¬ 
ous  frequencies.  The  interaction  between  an  imposed  monochromatic,  time-dependent  acoustic 
disturbance  and  a  stea.dy  mean  shear  flow  in  a  two-dimensional  duct  was  studied  using  direct  nu¬ 
merical  simulation.  Unlike  previously  reported  numerical  results,  the  long  time  acoustic  response 
is  captured  through  implementation  of  accurate  non-reflecting  boundary  conditions.  Below  a 
certain  cut-off  frequency,  the  acoustic  field  is  a  nearly  planar  traveling  wave  propagating  axially 
along  the  duct.  Above  this  cut-off  frequency,  oblique  waves  are  generated  due  to  both  acoustic 
refraction  and  cross-stream  dependent  source  oscillation,  leading  to  an  alternating  pattern  of 
higher  acoustic  pressures  at  the  wall  and  the  centerline,  downstream  of  the  disturbance  source. 
At  resonant  conditions,  the  growth  of  the  oblique  wave,  which  is  nearly  transverse,  dominates 
the  axial  wave,  leading  to  a  phase  chcinge  of  180  degrees  after  their  interaction.  The  thickness  of 
the  acoustic  boundary  layer  and  its  response  to  imposed  disturbances  are  also  in  good  agreement 
with  theory.  These  results  are  consistent  with  previously  published  theoretical  predictions,  and 
represent  their  first  numerical  verification.  The  code  is  being  modified  to  include  nonpremixed 
combustion,  which  will  allow  studying  acoustic  interactions  with  a  chemically  reactive  mean 
flow.  The  manuscript  in  the  Appendix  is  to  appear  in  the  AIAA  Journal  either  in  1995  or  1996, 
and  represents  most  of  the  work  done  in  connection  with  this  contract. 
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1.  INTRODUCTION  AND  RESEARCH  OBJECTIVES 


Combustion-driven  acoustic  instabilities  in  solid  rocket  motors  arise  as  a  result  of  a  complex 
interaction  between  the  acoustics  of  the  flow  field  and  oscillations  in  propellant  burning  rate. 
The  net  energy  associated  with  acoustic  instabilities  is  determined  by  a  balance  of  competing 
processes  that  feed  and  extract  energy  from  the  acoustic  field.  The  relative  importance  of  these 
processes  is  strongly  dependent  on  both  the  mode  of  oscillation  and  type  of  propellant,  thus 
making  it  difficult  to  isolate  physical  mechanisms  associated  with  individual  processes  when  the 
problem  is  studied  as  a  whole. 

The  focus  of  the  present  research  is  on  the  effects  of  acoustic  refraction  caused  by  the  in¬ 
teraction  of  disturbances  imposed  on  a  viscous  mean  flow  through  a  two  dimensional  channel. 
This  is  accomplished  through  Direct  Numerical  Simulation  (DNS)  of  the  flow  field  in  several 
model  problems  tailored  to  address  specific  features  of  the  interaction  between  the  mean  flow 
and  the  acoustic  disturbances.  The  overall  objective  is  to  develop  a  methodology  that  will  allow 
one  to  extract  acoustic  information  from  data  generated  through  DNS.  It  is  expected  that  this 
process  will  lead  to  a  better  understanding  of  acoustic-mean  flow  interactions  in  solid  rocket 
motors.  Our  specific  objective  is  to  validate  the  numerical  simulations  with  analytical  results 
obtained  by  Wang  and  Kassoy  (1991),  thereby  developing  and  fine  tuning  a  procedure  that  can 
be  used  to  track  energy  in  the  mean  flow  for  various  acoustic  modes  of  oscillation.  To  date, 
we  have  modified  an  existing  DNS  computer  code  and  have  obtained  results  for  unsteady  flow 
in  a  rectangular  duct.  The  DNS  procedure  is  described  in  Section  2.  Our  progress  to  date  is 
presented  in  Section  3.  Finally,  an  outline  of  our  future  research  plans  is  provided  in  Section  4. 

2.  DIRECT  NUMERICAL  SIMULATION  PROCEDURE 

A  direct  numerical  simulation  code  for  2D  compressible  flow  developed  by  Baum  (1993)  wa^ 
modified  and  is  used  in  our  investigations.  The  continuity,  momentum  and  energy  equations  for 
compressible  flow  are  solved  numerically.  Pressure,  density,  and  temperature  are  related  through 
the  ideal  gas  equation  of  state.  The  specific  heats  are  presumed  constant  and  fluid  properties 
such  as  dynamic  viscosity  and  thermal  conductivity  are  assumed  to  be  functions  of  temperature. 
The  equations  are  solved  using  a  sixth-order  accurate  compact  finite  difference  scheme  (Lele, 
1992)  for  evaluating  spatial  derivatives.  Time  advancement  is  achieved  through  a  third-order 
accurate  Runge-Kutta  scheme.  The  Navier  Stokes  Characteristic  Boundary  Condition  (NSCBC) 
procedure  developed  by  Poinsot  and  Lele  (1992)  is  used  to  implement  boundary  conditions  at 
the  lateral  and  streamwise  edges  of  the  computational  domain. 

Figure  1  is  a  schematic  representation  of  our  computational  domain  and  coordinate  system 
definition.  The  aspect  ratio  is  defined  as  L/2/i.  The  grid  spacing  in  the  computational  mesh 
is  varied  depending  on  the  computational  resolution  required,  and  to  evaluate  the  sensitivity  of 
the  analysis  to  grid  spacing.  The  equations  are  cast  in  dimensionless  form  and  time-advanced 
to  obtain  either  steady  or  unsteady  solutions.  The  duct  length  L,  the  speed  of  sound  in  the  far 
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field  Coo»  aJid  a  characteristic  acoustic  time,  based  on  L  and  Coo  sire  used  as  reference  length, 
velocity,  and  time  scales.  Other  variables  are  rendered  dimensionless  by  their  properties  in  the 
fax  field. 

3.  PROGRESS  TO  DATE 

As  reported  previously,  several  test  cases  were  run  to  validate  the  modifications  that  were 
made  to  the  2D  code,  including  the  behavior  of  the  nonrefiecting  boundary  conditions,  and  the 
overall  behavior  of  the  code  for  stability  and  convergence.  A  critical  requirement  for  accurately 
capturing  acoustic  wave  structure  is  a  nonreflecting  outlet  boundary,  where  waves  generated  at 
the  inlet  can  paiss  through  with  little  or  no  reflection  back  into  the  computational  domain. 

After  validating  the  boundary  condition  treatment  and  the  numerical  scheme,  perturbations 
were  imposed  on  the  inlet  axial  velocity.  In  order  to  compare  numerical  results  with  the  analytical 
results  of  Wang  and  Kassoy  (1992),  it  is  necessary  to  impose  a  plane  wave  (constant  amplitude) 
disturbance.  However  it  is  not  possible  to  implement  this  form  of  perturbation  in  the  numerical 
simulation  because  of  the  no  slip  condition  that  is  required  at  the  lateral  boundaries.  The  actual 
form  of  the  disturbance  used  in  the  simulation  is  shown  in  Figure  2,  and  is  mathematically 
expresed  as 


u  z=  AU{y)  sin  [Stt  (y  +  0.1)]  sin(fit),  (1) 

where  u  is  the  axial  velocity  disturbance,  ft  is  the  dimensionless  frequency  of  the  disturbance 
defined  through  ft  =  where  ft'  is  the  dimensional  frequency,  and  a  is  the  speed  of  sound. 

The  advantage  of  this  amplitude  structure  is  that  it  does  not  introduce  any  discontinuity  or 
sharp  gradients  so  that  relatively  large  values  for  A  may  be  used.  Also,  by  performing  a  Fast 
Fourier  Transform  (FFT)  in  y  for  a  discrete  representation  of  Equation  1,  it  is  found  that  the 
dominant  transverse  spatial  components  are  the  zeroth  and  first  modes,  with  higher  modes  being 
negligibly  small. 

With  an  expression  for  the  imposed  disturbance  of  the  axial  velocity  specified,  the  acoustic 
pressure  field  under  various  driving  frequencies  was  studied.  Several  values  of  angular  frequency 
were  considered.  Guided  by  the  work  of  Wang  and  Kassoy  (1992),  these  values  were  chosen  to 
represent  typical  characteristics  of  different  frequency  regimes.  In  all  cases,  the  computations 
were  performed  for  up  to  four  acoustic  time  units,  compared  to  only  one  acoustic  time  unit 
in  previous  numerical  simulations. 

1.  Acoustic  pressure  field  for  ft  =  10.  This  is  the  lowest  frequency  studied  in  the  simulations, 
and  corresponds  to  the  dimensionless  frequency  ft^y  =  2  used  in  the  analysis  of  Wang  and 
Kassoy  (1992).  Figures  3(a)  and  (b)  represent  the  total  acoustic  pressure  time  history  and 
pressure  contours  at  the  indicated  locations.  Notice  that  in  Figure  3(a)  the  amplitude  of 
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the  acoustic  pressure  at  the  centerline  is  slightly  smaller  than  that  at  the  wall.  This  is  a 
consequence  of  acoustic  refraction  that  deflects  the  acoustic  waves  toward  the  boundaries. 
The  constant  wave  amplitude  shown  in  Figure  3(a)  indicates  that  the  only  wave  existing 
in  the  duct  is  that  imposed  at  the  inlet.  The  contour  plot  in  Figure  3(b)  confirms  this  by 
showing  that  the  acoustic  pressure  gradient  is  primarily  in  the  axial  direction. 

2.  Acoustic  pressure  field  for  =  40.  The  acoustic  pressure  time  history  at  fixed  locations 
is  presented  in  Figure  4  for  a  dimensionless  forcing  frequency  H  =  40.  Compared  with 
Figure  3(a),  it  reveals  several  interesting  features.  Initially,  after  the  arrival  of  the  acoustic 
waves,  an  increase  in  their  amplitudes  is  observed.  The  amplitudes  then  staxt  to  oscillate 
over  time.  This  variation  in  amplitude  can  be  explained  by  the  analytical  work  of  Wang 
and  Kctssoy  (1992).  They  identified  a  cutoff  frequency  Ct  =  5n7r,  n  =  1, 2, 3, ...  below  which 
the  waves  remain  planar  in  the  duct.  For  a  plane  wave  with  above  IOtt,  oblique  waves 
axe  generated  as  the  disturbance  moves  downstream.  The  path  of  the  oblique  waves  can 
be  traced  by  examining  the  acoustic  pressure  contours  in  Figure  5.  The  direction  of  the 
maximum  pressure  gradient  is  an  indication  of  the  direction  the  waves  propagate.  The 
wave  front  is  aligned  normal  to  this  direction.  The  intersection  between  the  developing 
oblique  waves  and  the  axially  travelling  waves  results  in  the  variation  of  amplitude. 

Another  important  result  may  be  obtained  by  examining  the  acoustic  pressure  time  his¬ 
tory  along  the  wall  and  at  the  centerline  at  several  streamwise  stations.  These  data  are 
presented  in  Figure  6.  A  consequence  of  acoustic  refraction  due  to  shear  flow  interaction 
is  that  the  acoustic  pressure  is  higher  at  the  wall  than  at  the  centerline  for  downstream 
propagating  waves.  However,  in  the  presence  of  oblique  waves  this  is  not  true  everywhere 
in  the  domain.  Instead,  by  studying  Figures  5  and  6,  one  notes  that  higher  acoustic  pres¬ 
sure  exists  only  in  those  places  where  oblique  waves  pass  through.  This  further  indicates 
the  importance  of  oblique  waves  in  the  study  of  the  acoustic  field, 

3.  Acoustic  pressure  field  for  fi  =  IOtt.  This  angular  frequency  corresponds  to  the  second 
resonant  mode  of  the  channel.  Figures  7(a)-(h)  depicts  the  acoustic  pressure  time  his¬ 
tory  at  the  centerline  and  at  the  wall  at  different  axial  locations  along  the  channel.  The 
general  feature  of  the  resonant  oscillation  is  displayed  at  locations  (Ic)  and  (Iw),  where 
the  acoustic  pressure  amplitudes  increase  with  time.  Figures  7(b)  and  (c)  depict  slightly 
different  features.  One  can  observe  that  at  some  locations  the  signal  initially  decreases 
and  subsequently  increcises.  A  similar  behavior  is  presented  in  Wang  and  K2tssoy  (1992). 
More  interesting  behavior  can  be  seen  further  downstream.  After  the  initial  decrease  in 
amplitude  due  to  destructive  interference,  the  wave  amplitudes  at  all  locations  start  to 
increase  with  time.  Then  at  some  locations  the  wave  amplitudes  remain  nearly  constant, 
suggesting  that  acoustic  energy  due  to  resonant  oscillations  is  somehow  being  transferred 
out  of  the  acoustic  field  at  these  locations.  However,  at  other  locations  the  wave  ampli¬ 
tudes  start  to  decrease  and  eventually  reach  a  value  close  to  zero.  After  that,  the  pressure 
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amplitudes  increase  once  again  having  undergone  a  phase  change  of  tt  with  respect  to  the 
phase  at  the  inlet.  This  phase  change  is  illustrated  in  Figure  8,  where  data  from  locations 
(6c)  and  (6w)  have  been  overlaid.  A  similar  behavior  was  suggested  by  Wang  and  Kassoy 
(1992);  however,  since  previous  numerical  simulations  were  limited  to  a  few  wave  cycles 
only,  the  results  had  not  been  verified.  The  present  numerical  simulations  overcome  this 
limitation,  and  the  results  show  good  qualitative  agreement  with  analytical  predictions. 


Acoustic  Boundary  Layer  Analysis.  The  existence  of  an  acoustic  boundary  layer  is  another 
interesting  feature  that  arises  due  to  the  acoustic-mean  flow  interaction.  Represented  in  Figures 
9(a)-(e),  it  is  the  region  close  to  the  wall  within  which  an  overshoot  of  axial  acoustic  velocity 
can  be  observed  The  magnitude  of  the  overshoot  reaches  a  maximum  value  whenever  the  axial 
velocity  in  the  core  changes  sign.  This  effect  is  called  Richardson’s  annular  effect  (Baum  and 
Levine,  1978,  and  Wang  and  Kassoy,  1992).  The  thickness  of  the  acoustic  boundary  layer  is 
given  in  analytical  form  by 


6  =  5 


f  M 
\QReJ 


(2) 


where  M  is  the  maximum  centerline  Mach  number  and  Re  is  the  Reynolds  number  based  on 
the  centerline  velocity  and  duct  width.  For  the  field  shown  in  Figures  9(a)-(e)  with  Re  =  320 
and  M  w  0.8,  Equation  2  predicts  an  acoustic  boundary  layer  thickness  of  5.6 

Preliminary  Results  for  Flow  with  Sidewall  Injection.  For  this  simulation,  the  left  boundary 
of  the  computational  domain  is  closed;  i.e,  flow  is  not  permitted  through  the  boundary.  Flow 
through  the  channel  is  generated  by  specifying  mass  injection  from  the  lateral  boundaries  along 
y  =  ±h  (see  Figure  1).  The  boundary  conditions  imposed  are  (1)  u  =  t;  =  0,T  =  const,  along 
a;  =  0,  (2)  nonreflecting  boundary  at  outlet,  and  (3)  tt  =  0,t;  =  ±Vb,T  =  const,  along  y  =  ±/i, 
where  Vq  is  the  magnitude  of  the  injection  velocity,  taken  to  be  0.01  in  the  present  simulation. 
Figures  10(a)  and  (b)  show  the  comparison  between  the  numerical  and  inviscid  analytical  steady 
flow  model.  The  difference  between  the  two  are  probably  due  to  the  viscous  nature  of  the  fluid 
assumed  in  the  numerical  simulation. 

After  the  steady  state  is  obtained,  a  small  pressure  disturbance  is  imposed  at  the  outlet 
through 


Ptotal, outlet  —  Psteady, outlet  ^P 


(3) 


where  e  is  a  asymptotic  parameter  taking  the  value  0.01  and 


p  =  Asin(Q,t) 


(4) 
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where  A,  the  magnitude  of  the  disturbance,  is  taken  as  the  steady  state  pressure  at  the  center 
of  the  outlet.  Dimensionless  frequencies  assumed  in  the  analysis  were  Cl  =  5, 10, 107r,40.  The 
vortical  component  of  the  axial  velocity  is  obtained  through  the  following  relation  suggested  by 
Price  and  Flandro  (1992)  and  Zhao,  et.  al.  (1994): 


Cortical  —  '^^steady  l^center 


(5) 


Figure  11  represents  the  comparison  between  «vorticai  profiles  for  the  different  frequencies.  The 
acoustic  boundary  layer  is  extended  towards  the  core  region  of  flow  due  to  the  convective  effect 
of  mass  injection.  However,  the  thickness  of  this  layer  is  still  dependent  on  the  frequency  as  it 
is  in  the  case  for  flow  without  mass  injection. 

4.  FUTURE  WORK 

We  are  currently  enhancing  the  2D  code  in  order  to  simulate  a  2D  nonpremixed  flame 
confined  between  the  walls  of  a  narrow  horizontal  duct  with  open  ends.  Single  step,  finite  rate 
chemistry  is  assumed  for  the  combustion  process,  with  the  unmixed  reactants  diffusing  into  the 
flame  zone  from  opposite  sides  of  the  channel.  Eventually  one  end  of  the  channel  will  be  closed, 
and  anoustic  perturbations  will  be  imposed  on  the  reacting  flow  field.  This  is  the  first  step 
toward  modeling  the  effects  of  acoustics  on  the  combustion  response  of  a  nonpremixed  flame 
confined  near  the  surface  of  a  burning  propellant.  In  the  future,  this  analysis  will  be  extended 
to  simulate  the  response  of  a  premixed  flame  adjacent  to  the  burning  propellant. 
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Figure  1  Schematic  representation  of  computational  domain  and  coordinate  system  defi¬ 
nition. 


Figure  2  Profile  of  amplitudes  of  axial  velocity  disturbance  used,  normalized  with  respect 
to  the  centerline  velocity  for  A  —  1. 
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Total  acoustic  pressure 


Acoustic  time 


Figure  4  Total  acoustic  pressure  time  history  at  x  =  0.15,  center  (dashed)  and  wall  (solid) 
for  =  40,  mesh  81  X  81  and  acoustic  Re  =  20,000. 


y 


(b) 


Figure  10  (a)  Steady  state  axial  velocity  profiles  at  i  =  0.5  for  acoustic  Re  =  20,000 

with  side  wall  injection.  Solid  line:  analytical  solution,  Dashed  line:  numerical 


solution;  (b)  Steady  state  transverse  velocity  profile. 


axial  velocity 


Figure  11 


Profiles  of  the  vortical  component  of  the  axial  velocity  at  x  =  0.5  for  acoustic 

ile  =  20, 000  with  side  wall  injection.  - =  5,  -  = 

10,  - :fl  =  IOtt,  -  :Sl  =  40 
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Abstract 

The  interaction  between  an  imposed  monochromatic,  time-dependent  acoustic  disturbance 
and  a  steady  mean  shear  flow  in  a  two-dimensional  duct  is  studied  using  direct  numerical 
simulation.  Unlike  previously  reported  numerical  results,  the  long  time  acoustic  response 
is  captured  through  implementation  of  accurate  non-reflecting  boundary  conditions.  Below 
a  certain  cut-off  frequency,  the  acoustic  field  is  a  nearly  planar  traveling  wave  propagating 
axially  along  the  duct.  Above  this  cut-off  frequency,  oblique  waves  are  generated  due  to  both 
acoustic  refraction  and  cross— stream  dependent  source  oscillation,  leading  to  an  alternating 
pattern  of  higher  acoustic  pressures  at  the  wall  and  the  centerline,  downstream  of  the  dis¬ 
turbance  source.  At  resonant  conditions,  the  growth  of  the  oblique  wave,  which  is  nearly 
transverse,  dominates  the  axial  wave,  leading  to  a  phase  change  of  180  degrees  after  their 
interaction.  The  thickness  of  the  acoustic  boundary  layer  and  its  response  to  imposed  distur¬ 
bances  are  also  in  good  agreement  with  theory.  These  results  are  consistent  with  previously 
published  theoretical  predictions,  and  represent  their  first  numerical  verification. 
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Nomenclature 


ttoo  reference  speed  of  sound,  also  reference  velocity 
Cp  constant  pressure  specific  heat,  temperature  independent 
ct  internal  energy  per  unit  volume 

h  duct  half-width 

iyjyk  denote  coordinate  directions  in  Cartesian  tensor  notation 
L  computational  length  of  duct 

M  mean  flow  Mach  number 

p  pressure 

Pr  Prandtl  number 

qj  heat  flux  vector  component  in  the  j— th  direction 
Re  acoustic  Reynolds  number  in  the  calculation 
Rcc  duct  Reynolds  number 

T  temperature 

t  time  variable 

Uk  velocity  component  in  the  fc-th  coordinate  direction 
u  axial  acoustic  velocity  disturbance 
V  velocity  component  in  the  y  coordinate  direction 
Xyy  spatial  variables  along  the  x  and  y  coordinate  directions 
Xj  independent  spatial  variable  along  the  j-th  coordinate  direction 

Greek  symbols 

a  computational  duct  aspect  ratio 

7  constant  specific  heat  ratio 

8  acoustic  boundary  layer  thickness 

€  small  parameter  used  in  perturbation  analysis 
K  viscosity  power  law  exponent 
coefficient  of  dynamic  viscosity 
u  coefficient  of  kinematic  viscosity 

fl  dimensionless  frequency  of  imposed  oscillation 
Tij  stress  tensor 

Subscripts  and  superscripts 

’  denotes  dimensional  quantity 

oo  denotes  reference  quantity,  chosen  at  stagnant  conditions 
denotes  acoustic  disturbance 
s  undisturbed  steady  state 

Introduction 

The  interaction  of  an  acoustic  wave  disturbance  with  a  shear  flow  provides  a  mechanism  for 
transfer  of  energy  between  the  mean,  and  various  modes  of  the  acoustic  flow.  This  problem 
is  significant  sls  it  provides  a  vital  link  in  the  chain  of  events  that  could  lead  to  combustion- 
driven  acoustic  instability  in  solid  rocket  motors.  However,  in  practical  motors,  the  various 
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mechanisms  that  contribute  to  amplification  and  damping  of  acoustic  energy  are  intimately 
related  to  the  acoustic  mode,  and  the  type  of  propellant  lining  the  motor  surface.^  This  makes 
it  difficult  to  isolate  physical  mechanisms  associated  with  individual  processes  when  the  full 
problem  is  investigated.  Thus,  the  focus  of  this  paper  is  on  the  mechanism  of  energy  exchange 
between  the  mean  and  various  modes  of  the  acoustic  flow  in  a  sufficiently  simplified  flow 
situation  for  which  theoretical  results  are  available.  This  problem  is  investigated  via  Direct 
Numerical  Simulation  (henceforth,  DNS)  of  the  interaction  between  an  imposed  acoustic 
velocity  disturbance  in  an  otherwise  steady  shear  flow  established  in  a  two-dimensional  duct. 
In  DNS,  the  complete  time-dependent  system  of  equations  is  solved  without  any  explicit 
modeling,  such  as  turbulence  modeling.  This  procedure  ha>s  been  successfully  applied  to  a 
variety  of  incompressible,  compressible,  and  reacting  flows. ^ 

Much  of  past  analytical  work  has  focused  on  obtaining  quasi-steady  solutions  to  a  lin¬ 
earized  wave  equation  describing  acoustic  wave  propagation  in  a  fully  developed  duct  flow.^ 
Using  this  as  a  basis,  both  downstream  propagating^  and  upstream  propagating  acoustic 
waves'^*®  have  been  investigated.  Results  suggest  that  in  the  former  case,  the  acoustic  pres¬ 
sure  at  the  wall  is  larger  than  that  at  the  centerline;  the  reverse  is  true  in  the  latter  case. 
Since  this  method  presumes  the  form  of  the  solution,  it  cannot  predict  the  evolution  from 
an  initial  disturbance  to  the  quasi-steady  form.  Baum  and  Levine^,  and  subsequently,  Wang 
and  Kassoy^  proceeded  to  rectify  this  shortcoming  by  investigating  an  initial-boundary  value 
problem  albeit  using  entirely  different  approaches.  Baum  and  Levine  solved  the  compressible 
Navier-Stokes  equations  numerically  for  a  Reynolds  number,  based  on  centerline  velocity  and 
duct  width,  of  approximately  10^.  A  turbulence  model  was  included  in  their  calculations.  A 
disturbance  was  introduced  at  a  certain  cross  section  and  its  evolution  was  tracked.  They 
found  that  acoustic  refraction  effects  increased  with  the  frequency  of  the  imposed  distur¬ 
bance.  The  relationship  between  the  wall  and  centerline  acoustic  pressure  for  upstream  and 
downstream  propagation  of  acoustic  waves,  in  a  qualitative  sense,  was  consistent  with  that 
predicted  by  the  quasi-steady  theory. At  the  outflow  boundary,  they  had  difficulty  in  pre¬ 
scribing  non-reflecting  conditions,  and  hence  they  terminated  their  computations  when  the 
imposed  wave  reached  the  downstream  boundary  of  their  computational  domain.  Thus  the 
long-time  transient  solution,  which  has  since  been  demonstrated  to  be  extremely  important'^ 
was  not  available  in  their  work.  This  may  also  be  the  reason  why  they  did  not  observe  oblique 
waves  and  resonant  mode  oscillations.  Using  a  rational  perturbation  procedure,  Wang  and 
Kassoy*^  solved  an  initial-boundary  value  problem  to  describe  acoustic  processes  in  a  two- 
dimensional  shear  flow  in  a  duct  for  small  mean  flow  Mach  number  and  large  Reynolds 
number.  Their  results  indicate  that  as  a  result  of  leading  order  axial  wave  refraction  by 
the  shear  flow,  oblique  propagating  waves  develop.  This  refraction  effect  increases  with  the 
driving  frequency  of  the  imposed  acoustic  disturbance  and  the  mean  flow  Mach  number.  The 
oblique  waves  evolve  into  purely  amplifying  transverse  waves  at  frequencies  corresponding  to 
resonance.  The  initial-boundary  value  approach  discussed  above  differs  from  what  may  be 
termed  a  normal  mode  analysis  that  is  often  applied  to  problems  involving  combustion-driven 
acoustic  instability  in  solid  rocket  motors®»^»^®  and  pulsed  combustors. In  the  traditional 
approach, the  amplitudes  of  either  the  classical  acoustic  modes  or  modes  suitably  mod- 
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ified  to  account  for  processes  that  influence  the  bztsic  acoustic  modes  are  described  via  an 
infinite  system  of  coupled  ordinary  differential  equations.  The  idea  is  then  to  obtain  con¬ 
ditions  that  could  lead  to  growth  or  amplification  of  different  modes,  and  also  to  predict 
conditions  that  could  lead  to  limit-cycle  type  behavior. 

The  objective  in  the  present  paper  is  to  verify  the  existence  of,  and  to  understand  the 
role  of  oblique  waves  generated  through  acoustic  refraction  when  a  monochromatic,  acoustic 
velocity  disturbance  introduced  at  a  fixed  duct  location  is  allowed  to  interact  with  a  steady 
shear  flow  in  a  two-dimensional  duct.  In  particular,  we  seek  to  validate  the  analytical  re¬ 
sults  obtained  by  Wang  and  Kassoy^  and  use  their  predictions  to  examine  and  interpret  the 
complex  acoustic  structure  that  is  obtained.  The  approach  adopted  is  an  accurate  solution 
to  the  compressible  Navier-Stokes  equations  through  DNS  of  the  problem. 

Problem  Formulation 

The  full  compressible  two-dimensional  Navier-Stokes  equations  are  solved  as  an  initial  value 
problem  using  a  code  originally  developed  by  Baum  and  Poinsot^^  for  reacting  flows  and 
adapted  for  the  present  research.  The  computational  domain  and  coordinate  system  are 
represented  in  Fig.  1.  The  equations  are  rendered  dimensionless  using  the  duct  length  L, 
speed  of  sound  at  stagnant  conditions  aoo>  the  corresponding  density  /9oo>  and  the  dynamic 
pressure  Poo®^)-  The  reference  temperature  is  (7  —  l)Too  The  reference  viscosity  used  is 
Hoo  =  ^i(Too)^  The  reference  time  is  thus  the  acoustic  time  I/Ooo.  The  dimensionless 
parameters  governing  the  problem  are  the  Reynolds  number  defined  as 


Re= 

fJ'OO 


(1) 


and  the  Prandtl  number  Pr  =  fiooCpfXooi  where  Ao©  =  A  (To©).  Note  that  the  duct  aspect 
ratio  defined  as 


(2) 


although  not  a  paxameter,  appears  in  the  analysis  since  the  computational  domain  is  finite. 
By  applying  non-~reflecting  boundary  conditions  as  discussed  in  the  next  section,  the  flow  field 
in  a  duct  of  infinite  length  is  effectively  simulated.  Henceforth  all  quantities  are  dimensionless. 
Dimensional  quantities  where  appropriate  are  indicated  by  a  prime.  The  dimensionless  form 
of  the  governing  equations  written  in  Cartesian  tensor  notation  appears  below: 

Continuity, 
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momentum, 


dpui 

dt 


dpUiUj 

dxj 


dp  1  dTjj 
dxi  Re  dxj  ’ 


and  energy, 

dpet  djpet  +  p)uj  _  1  d^UjTjj)  _  1  dqj 

dt  dxj  ~  Re  dxj  RePrdxj' 


(4) 


(5) 


The  internal  energy  per  unit  volume  is  given  by 


1  p 

=  ::pUkUk  H - 

z  If  ^ 


(6) 


The  stress  tensor  Tij  and  the  heat  flux  vector  qj  are  taken  for  a  Newtonian  fluid  with  the 
Fourier  model  for  heat  conduction  as  follows: 


Tij  =  p 


dui  duj  2^  duk 
dxj  ^  dxi  3  dxk^ 


Qi  =  -A 


dT 

dxi 


(7) 


The  dynamic  viscosity  p  is  dependent  on  temperature  such  that  p  ~  T'‘,  where  the  exponent 
K  is  a  constant.  The  thermal  conductivity  A  is  related  to  the  dynamic  viscosity  through  the 
constant  Prandtl  number.  The  equation  of  state  is 


p  =  pt'^!—^.  (8) 
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Numerical  Method  and  Boundary  Treatment 

A  compact  sixth-order  accurate  finite  difference  scheme  developed  by  Lele^^  is  used  to  ap¬ 
proximate  the  first  and  second  derivatives  in  the  governing  equations.  The  order  of  accuracy 
is  fourth  and  third  respectively  at  points  adjacent  to,  and  on  the  boundary  of  the  compu¬ 
tational  domain.  This  scheme  requires  solution  of  a  tridiagonal  matrix  in  order  to  compute 
the  derivatives.  The  advantage  of  the  compact  scheme  over  more  familiar  finite-difference 
schemes  is  that  it  provides  spectral-like  resolution  of  both  the  amplitude  and  phase  of  the 
solution,  thereby  enabling  the  scheme  to  accurately  capture  acoustic  wave  behavior.  A  com¬ 
plete  analysis  of  this  scheme  including  a  comparison  of  dispersion  errors  with  traditional 
finite  difference  schemes  is  presented  by  Lele.^^  A  third-order  explicit  Runge-Kutta  scheme 
is  adopted  for  time-advancement  of  the  semi-discrete  equations. 
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One  critical  requirement  for  capturing  the  wave  characteristics  of  the  acoustic  field  is 
accurate  treatment  of  boundary  conditions.  In  the  present  paper,  a  systematic  boundary 
treatment  method,  namely  Navier-Stokes  Characteristic  Boundary  Condition  (NSCBC),  is 
utilized.^^  Based  on  the  idea  of  characteristics  for  Euler  equations^®  the  NSCBC  assumes  that 
the  waves  propagated  by  Navier-Stokes  equations  are  associated  with  only  the  hyperbolic 
part  of  the  equation. Thus  the  first  derivative  normal  to  a  boundary  appearing  in  the 
convective  term  of  the  Navier-Stokes  equations  may  be  rewritten  in  terms  of  the  amplitude 
variations  of  the  characteristic  waves  represented  by  Navier-Stokes  equations.  For  those 
characteristic  waves  propagating  out  of  the  domain,  the  values  of  their  amplitude  variation  are 
completely  defined  by  data  from  within  the  computational  domain  through  use  of  one-sided 
approximations.  Numerical  stability  is  achieved  since  this  amounts  to  upwind  differencing. 
For  those  waves  propagating  into  the  computational  domain,  boundary  conditions  are  needed 
to  specify  the  values  of  their  amplitude  variations.  There  is  no  exact  method  to  compute 
the  amplitude  variations  of  these  incoming  waves.  A  set  of  locally  one-dimensional  inviscid 
(LODI)  equations  are  used  to  infer  the  values  of  the  incoming  amplitude  variations  from  the 
boundary  conditions.  Derivatives  parallel  to  the  boundary  pose  no  problems  and  are  treated 
exactly  as  in  the  interior.  The  treatment  of  viscous  conditions  is  based  on  the  method 
suggested  in  Ref.  14.  These  terms  go  to  zero  as  the  viscosity  goes  to  zero.  For  a  detailed 
discussion  of  the  NSCBC  procedure,  the  paper  by  Poinsot  and  Lele^^  may  be  consulted. 

In  this  paper  the  flow  considered  is  subsonic  with  a  mean  flow  Mach  number  M  «  0.08. 
Unless  otherwise  noted,  a  =  5^  Re  =  20000,  and  Pr  =  1.  Thus  the  Reynolds  number  baised 
on  the  mean  flow  velocity  and  duct  width,  Rcc  =  MRefa  =  320. 

Reflecting  boundary  conditions  are  used  at  the  inlet.  The  transverse  velocity  compo¬ 
nent  V  is  set  to  zero,  the  axial  velocity  profile  u  =  u(0,y,t),  such  that  u(0,— l/2a,t)  = 
ti(0,-f-l/2a,t)  =  0,  is  prescribed,  and  a  constant  inlet  temperature,  T(0,j/,t)  =  1/(7  —  1) 
is  prescribed.  At  the  lateral  boundary,  no-slip  velocity  and  isothermal  conditions  are  pre¬ 
scribed.  At  the  outflow,  non-reflecting  boundary  conditions  baised  on  the  method  developed 
in  Ref.  14  are  imposed.  This  allows  acoustic  waves  to  propagate  out  of  the  domain  with  little 
or  no  reflections,  and  thus  enables  us  to  carry  out  long-time  stable  computations.  A  steady 
flow  is  established  by  initializing  the  domain  with  stagnant  flow,  and  raising  the  pressure  at 
the  inlet  by  5%  of  the  initial  uniform  pressure,  p  =  1,  This  steady  state  flow  field  is  used  as 
initial  conditions  for  several  transient  problems  in  which  a  time-dependent  perturbation  is 
imposed  on  the  incoming  streamwise  velocity. 

Results  and  Discussion 

Code  Validation 

Several  test  problems  were  simulated  to  verify  that  the  code  and  boundary  condition  treat¬ 
ment  are  satisfactory.  Both  non-reflecting  and  reflecting  boundary  conditions  were  examined.^® 
In  order  to  test  the  NSCBC  procedure,  a  very  viscous  test  problem^^  was  simulated.  The 
inlet  conditions  imposed  are 


u(0,  y,  t)  -  uo  cos^  (rya) ,  t)(0,  y,  t)  =  0,  T(0,  y,  t)  =  To,  (9) 
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where  uq,  Tq  are  constants.  Non-reflecting  conditions  are  imposed  at  the  outlet.  The 
Reynolds  number  Re  =  2000  and  a  uniform  41  X  41  mesh  was  used.  Steady  state  con¬ 
tour  plots  of  u/uo,  100(p  -  Poo)/Poo,  and  100(r  -  Too)IT^  are  virtually  identical  to  past 
numerical^^  and  analytical  results,^"*  thereby  validating  the  code  and  boundary  conditions. 
For  the  sake  of  brevity,  in  Fig.  2,  only  contour  plots  of  temperature  are  presented.  Tran¬ 
sient  results  presented  in  the  remainder  of  this  section  were  obtained  on  a  uniform  81  X  81 
mesh.  Grid  independence  was  verified  by  repeating  the  calculations  on  a  161  X  161  grid.  AU 
calculations  were  performed  on  a  Cray  C-90  at  the  San  Diego  Supercomputing  Center. 

Nature  of  imposed  disturbance 

A  transient  acoustic  field  is  generated  by  imposing  a  single  frequency  disturbance  on  the 
inlet  axial  velocity.  No  slip  at  the  lateral  walls  demands  that  the  amplitude  of  the  disturbance 
be  zero  at  p  =  ±l/2a.  Although  the  analysis  in  Ref.  7  is  quite  general,  the  results  they  present 
are  for  constant  amplitude  and  linear  variation  of  amplitude  and  thus  our  results  can  only 
be  qualitatively  related  to  their  theoretical  results.  The  axial  velocity  disturbance  u{y,t)  was 
chosen  to  have  the  following  mathematical  form 


u{y,t)  =  A{7(y)sin 


sin  fit 


(10) 


where  fi  is  the  dimensionless  frequency  defined  through  fi  =  Cl'Lfaoo,  and  U{y)  is  the  initial, 
steady  axial  velocity  at  the  inlet.  The  quantity  A  is  a  constant  that  controls  the  magnitude 
of  the  disturbance.  The  disturbance  is  introduced  by  specifying  an  axial  inlet  velocity  profile 
of  the  form 


w(0,P,t)  =  U(y)  +  €u(y,t) 


(11) 


The  quantity  e  is  the  small  parameter  introduced  in  Ref.  7.  By  performing  a  Fast  Fourier 
Transform  (FFT)  in  y  on  the  discrete  version  of  Eq.  (10),  it  is  found  that  the  primary 
transverse  mode  components  are  the  zeroth  and  first  modes,  with  higher  mode  amplitudes 
being  less  than  2%  of  these.  Thus  it  is  reasonable  to  expect  that  the  simulation  results  will 
compare  well  with  Ref.  7,  at  least  qualitatively. 

Theory^  for  p— independent  amplitude  function  suggests  that  for  Q,  <  amr,  n  =  1, 2, 3,-  • 
only  axial  waves  would  exist.  When  fl  >  anir,  refraction-induced  oblique  waves  begin  to 
appear.  In  the  present  simulation  a  =  5,  and  since  the  meetn  flow  is  symmetric  with  respect 
to  the  duct  center  plane,  only  even  values  of  n  are  relevant.  This  leads  to  the  first  cut-off 
frequency  of  IOtt  (corresponding  to  n  =  2)  above  which  oblique  waves  are  expected  to  be 
generated.  The  leading  order  solution  is  uninfluenced  by  the  mean  flow  except  through  a  bulk 
convection  effect.  Through  an  examination  of  the  second-order  acoustic  solution  Wang  and 
Kassoy^  predicted  that  fi  =  an-K  corresponds  to  resonant  conditions.  Classical  quasi-steady 
theory^  cannot  describe  resonant  oscillations,  and  previous  numerical  simulations  were  unable 
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to  capture  this  phenomenon  due  to  difficulties  in  implementing  outflow  boundary  conditions.® 
The  numerical  solutions  were  curtailed  when  the  disturbance  reached  the  outflow  boundary, 
and  hence  the  long-time  behavior,  required  to  observe  resonant  eflTects,  was  not  studied.®’^ 
In  the  present  work,  several  frequencies  were  examined.  In  the  following  sub-sections,  results 
are  presented  for  different  values  of  fi  that  were  chosen  to  represent  typical  characteristics  of 
different  frequency  regimes.  Since  the  no-slip  wall  condition  in  the  simulations  demands  that 
the  perturbation  amplitude  be  y-dependent,  it  is  important  to  recognize  that  higher  modes 
(non-planar)  are  expected  to  be  generated  directly  by  the  source,  in  addition  to  refraction  of 
the  planar  component.^  Using  the  decomposition  in  Ref.  7,  the  total  acoustic  pressure  p  and 
axial  acoustic  velocity  u  were  obtained  as 

p= -^\pix,y,t)-p,{x,y)],  u=^[u{x,y,t)- Us{x,y)]  (12) 

where  and  Us  are  the  undisturbed  steady  pressure  and  axial  velocity  respectively.  In  all 
crises,  computations  were  carried  out  up  to  four  acoustic  time  units,  with  e  =  0.1. 

Acoustic  Pressure  Field  for  fi  =  10 

The  frequency  =  10  is  below  the  first  “cut-ofP’  frequency  of  IOtt.  Fig.  3  represents  the 
total  acoustic  pressure  time  history  at  a:  =  0.6.  For  convenience,  the  time  scale  used  by  Wang 
and  Kassoy*^ 


(13) 


where  is  the  dimensional  time  is  also  shown  in  Fig.  3.  The  amplitude  of  acoustic  pressure  at 
the  centerline  is  slightly  smaller  than  that  at  the  wall  (also  true  at  other  x  locations)  due  to 
refraction  of  the  downstream  propagating  acoustic  wave.  An  examination  of  acoustic  pressure 
contours  suggests  that  sufficiently  far  downstream  from  the  inlet,  the  waves  are  nearly  planar. 
Thus  no  oblique  waves  appear  for  the  conditions  of  the  simulation.  This  result  is  in  qualitative 
agreement  with  analytical  prediction.^  An  FFT  analysis  of  the  acoustic  pressure  with  respect 
to  y  at  various  axial  locations  and  times  (see  Fig.  4)  shows  a  dramatic  decay  in  amplitudes 
of  the  first  transversal  spatial  component  while  very  little  decay  in  amplitudes  is  seen  for  the 
zeroth  component.  Since  the  latter  corresponds  to  the  planar  part  of  the  disturbance,  this 
result  is  consistent  with  the  earlier  remark  that  at  locations  sufficiently  far  downstream  from 
the  disturbance  source,  only  the  planar  component  survives. 

Acoustic  Pressure  Field  for  0  =  40 

The  forcing  frequency  fi  =  40  is  just  above  the  first  cut-off  frequency.  The  acoustic  pressure 
time  history  at  fixed  locations  is  presented  in  Fig.  5.  A  comparison  of  Fig.  5  with  Fig.  3  reveals 
several  interesting  features.  At  early  times,  after  the  first  arrival  of  the  waves,  an  increase 
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in  amplitude  can  be  observed.  Then  the  amplitudes  themselves  begin  to  oscillate  although 
at  a  much  reduced  frequency.  Recall  from  Ref.  7  that  a  plane  wave  (zeroth  transversal 
spatial  mode)  with  its  frequency  above  the  cut-off  frequency  generates  oblique  waves  as 
it  travels  downstream.  In  the  present  work,  the  acoustic  source  is  non-planar  due  to  the 
y— dependent  amplitude  function,  and  thus  oblique  waves  are  expected  to  be  excited,  in 
addition  to  refraction  of  the  axial  mode  giving  rise  to  oblique  waves.  The  result  of  these  two 
processes  is  a  variation  in  amplitudes  of  acoustic  pressure  as  depicted  in  Fig.  5.  Further  it 
is  evident  from  Fig.  5  that  the  pressure  amplitude  for  fi  =  40  is  greatly  enhanced  compared 
to  that  for  ft  =  10  (see  Fig.  3).  This  enhancement  arises  from  the  increased  amplitudes  of 
the  first  transversal  spatial  mode  as  can  be  seen  from  Figs.  6(a-b).  From  these  plots  it  is 
apparent  that  the  amplitudes  of  the  first  transversal  spatial  mode  grow  to  values  having  the 
same  order  as  those  of  the  zeroth  mode  while  the  amplitudes  of  the  zeroth  mode  varies  very 
little  compared  to  those  for  ft  =  10. 

Acoustic  Pressure  Field  for  ft  =  IOtt 

The  angular  frequency  ft  =  IOtt  corresponds  to  the  second  resonant  mode.  Figure  7  depicts 
the  acoustic  pressure  time  history  at  the  wall  and  centerline  at  different  axial  locations. 
The  general  feature  of  resonant  oscillation  is  displayed  at  locations  (Ic)  and  (Iw)  (here 
w  and  c  denote  wall  and  centerline  respectively),  where  the  acoustic  pressure  amplitudes 
increases  consistently  with  time.  However,  axial  stations  2  and  3  display  slightly  different 
features.  One  can  see  that  at  some  axial  locations  (stations  2w,  3w,  and  4c)  the  signal  initially 
decreases  and  subsequently  incre<Lses.  This  behavior  is  very  similar  to  that  reported  by  Wang 
and  Kassoy.^  Although  Wang  and  Kassoy^  focus  their  discussion  on  second  order  acoustic 
pressure,  their  qualitative  description  is  also  valid  for  the  total  anoustic  pressure  because  the 
first  order  acoustic  pressure  is  just  a  plane  sine  wave  propagating  axially,  for  y— independent 
source  amplitude.  The  acoustic  pressure  amplitude  at  different  axial  stations  increases  with 
time  after  sufficient  time  elapses.  However,  the  amplitude  of  acoustic  pressure  may  initially 
decrease  with  time  “owing  to  destructive  interference”.^  Previous  numerical  investigation  by 
Baum  and  Levine®  was  limited  to  very  short  time  solution  due  to  difficulties  alluded  to  earlier. 
Both  the  present  numerical  and  previous  analytical  results  show  that  an  examination  of  the 
short-time  solution  is  insufficient  to  conclude  whether  resonant  amplifications  axe  occurring. 

Acoustic  pressure  data  at  downstream  locations  (stations  4  through  8  in  Fig.  7)  reveals 
some  additional  features  that  have  not  been  reported  previously.  After  the  initial  decrease  in 
amplitude  due  to  “destructive  interference”  the  amplitudes  at  all  locations  start  to  increase 
with  time.  Then  at  some  locations  wave  amplitudes  quickly  stop  amplifying  and  remain 
neaxly  constant,  suggesting  that  extraction  of  acoustic  energy  (to  the  mean),  and  viscous 
dissipation  are  occurring.  On  the  other  hand,  at  other  locations  amplitudes  start  to  decrease 
again  and  eventually  reach  a  value  close  to  zero.  When  the  acoustic  pressure  amplitudes 
increase  once  again,  they  have  undergone  a  phase  change  of  tt  with  respect  to  the  phase  of 
the  signal  at  early  time.  This  can  be  seen  by  overlaying  the  data  from  locations  (6c)  and 
(6w)  of  Fig.  7,  not  included  here  for  the  sake  of  brevity.  This  behavior  of  acoustic  pressure 
waves  suggests  an  interaction  between  an  axially  traveling  wave  and  another  much  stronger 
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wave  which  arrives  after  the  axial  wave,  and  hajs  an  opposite  phase  with  respect  to  the  axial 
wave,  flawed  on  theoretical  results’^  it  may  be  concluded  that  the  alternating  appearance 
of  this  behavior  under  resonant  conditions  is  caused  by  the  interaction  of  axial  with  neaxly 
trajisverse  waves.  Unlike  the  ft  =  40  case,  transverse  waves  for  ft  =  IOtt  are  excited  and 
grow  with  time  due  to  resonance.  A  pictorial  representation  of  the  instantaneous  acoustic 
pressure  fields  is  presented  in  Figs.  8(a-d).  The  direction  along  which  the  pressure  gradient  is 
a  maximum  is  indicative  of  the  instantaneous  wave  structure  as  a  result  of  axial  and  transverse 
mode  interaction.  B£ised  on  a  comparison  of  the  present  results  with  theory*^,  the  following 
observations  may  be  made: 

1.  Present  simulation  shows  that  under  resonant  conditions,  significant  interciction  be¬ 
tween  the  nearly  transverse  and  axial  waves  only  appears  downstream,  sufficiently  far 
away  from  the  disturbance  source  probably  because  the  waves  need  a  certain  amount 
of  travel  distance  and  time  to  collect  sufficient  acoustic  energy  and  settle  into  nearly 
transverse  waves. 

2.  Prom  the  analysis  presented  in  Ref.  7,  where  the  source  amplitude  function  is  y— dependent, 
it  is  evident  that  the  magnitude  of  the  higher  modes  is  of  the  same  order  cis  the  axial 
mode  and  the  effect  of  refraction  of  the  axial  mode  is  of  lower  order.  This  suggests  that 
the  present  results  are  dominated  by  the  source  structure. 

Acoustic  Boundary  Layer 

The  existence  of  an  acoustic  boundary  layer  is  another  interesting  feature  that  arises  naturally 
due  to  the  acoustic-mean  flow  interaction.  It  is  a  region  close  to  the  wall  within  which  an 
overshoot  of  axial  acoustic  velocity  can  be  observed  a.s  depicted  in  Fig.  9.  The  magnitude 
of  the  overshoot  reaches  a  maximum  value  whenever  the  axial  acoustic  velocity  in  the  core 
changes  sign.  This  phenomenon  is  attributable  to  the  Richardson’s  annular  effect.®*^  The 
thickness  of  the  acoustic  boundary  layer  (J,  is  given  in  Ref.  7  by 
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This  thickness,  expressed  as  a  fraction  of  the  duct  width,  is  defined  as  the  distance  from 
the  wall  where  the  axial  velocity  reaches  97.3%  of  the  core  value.  Approximately  eight  grid 
points  are  present  within  this  boundary  layer,  a  value  that  is  sufficient  to  resolve  the  ax:oustic 
boundary  layer.  In  the  present  case,  6  =  5.6%  is  the  predicted  value.  Numerical  solution 
in  Figs.  9(a-d)  gives  a  thickness  of  5.0  to  5.5%  which  is  good  agreement  with  analytical 
prediction. 

Summary  and  Conclusions 

In  the  present  paper,  direct  numerical  simulation  was  performed  to  study  the  effect  of  an 
acoustic  disturbance  source  in  the  presence  of  a  mean  shear  flow  in  a  two-dimensional  duct. 
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Results  show  good  agreement  with  analytical  predictions  of  Wang  and  KassoyJ  Several  in¬ 
teresting  features  associated  with  this  problem  have  been  observed. 

1.  There  exists  a  critical  frequency,  above  which  oblique  waves  aje  generated.  The  interac¬ 
tion  between  oblique  waves  and  axially  traveling  waves  results  in  variation  of  amplitudes 
of  acoustic  pressure. 

2.  Amplified  oblique  waves  that  axe  nearly  transverse  due  to  resonant  disturbances  are 
observed  for  the  first  time.  The  rate  of  amplification  of  these  waves  is  much  higher 
than  that  of  the  axially  traveling  waves.  At  certain  locations  of  the  domain,  this  strong 
nearly  transverse  wave  interacts  with  the  axial  wave,  causing  a  dramatic  change  in  both 
wave  amplitude  and  pha.se. 

3.  For  resonant  disturbances  with  several  transversal  spatial  modes,  the  growth  rate  varies 
significantly  for  different  modes,  being  the  highest  for  the  first  transverse  mode. 

4.  The  present  simulation  resolved  the  acoustic  boundary  layer  very  well.  The  thickness 
estimated  through  numerical  simulation  is  in  agreement  with  that  predicted  by  theory. 
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Figure  Captions 


FIG.  1.  Schematic  representation  of  model  problem.  The  mean  flow  is  along  the  x'- 
axis.  Acoustic  perturbations  are  imposed  at  x'  =  0.  Non-reflectiong  boundary 
conditions  are  imposed  at  x'  =  L. 

FIG.  2.  Equally  spaced  contour  lines  of  100(!r  —  Too)lToo  at  t  =  25,  Re  =  2000 

FIG.  3.  Acoustic  pressure  time  history  at  x  =  0.6,  at  wall  (solid  line)  and  center 

(dashed  line),  fl  =  10,  mesh  =  81  X  81  and  Re  =  20000. 

FIG.  4.  Amplitudes  of  (a)  zeroth  and  (b)  first  transversal  spatial  mode  component  of 
acoustic  pressure  for  fi  =  10,  mesh  =  81  X  81  and  Re  =  20000. 

FIG.  5.  Acoustic  pressure  time  history  at  x  =  0.15,  center  (dashed)  and  wall  (solid) 
for  Q,  =  40,  and  Re  =  20000. 

FIG.  6.  Amplitudes  of  (a)  zeroth  and  (b)  first  transversal  spatial  mode  component  of 
acoustic  pressure  for  fl  =  40,  and  Re  =  20000. 

FIG.  7.  Acoustic  pressure  time  history  at  different  locations,  “c”:  centerline;  “w”: 

wall.  (1)  X  =  0.1,  (2)  X  =  0.2,  (3)  x  =  0.3,  (4)  x  =  0.4,  (5)  x  =  0.5,  (6) 
X  =  0.6,  (7)  X  =  0.7,  (8)  X  =  0.8  for  fl  =  IOtt,  and  Re  =  20000.  Note  acoustic 
time  scale  is  0  to  4  (0  to  140  in  the  time  scale  used  by  Wang  and  Kassoy), 
and  scale  for  acoustic  pressure  is  -5  to  5  in  all  the  figures 

FIG.  8.  Acoustic  pressure  contours  (equispaced)  at  (a)  t  =  3.0,  (b)  t  =  3.05,  (c) 
t  =  3.10,  (d)  t  =  3.14  for  =  IOtt,  and  Re  =  20000. 

FIG.  9.  Axial  acoustic  velocity  profiles  at  x  =  0.5  and  (a)  t  =  3.00,  (b)  t  =  3.14, 
(c)  t  =  3.30,  (d)  t  =  3.45,  (e)  t  =  3.62  for  fl  =  10,  mesh  =  81  x  161  and 
Re  =  20000.  Note  u  axis  scale  is  -0.6  to  0.6,  and  y  axis  scale  is  -0.1  to  -0.080. 
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